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I. INTRODUCTION 



The lattice approach to QCD pioneered by Wilson JJ and first realized numerically by 
Creutz 2! is based on the QCD action. Moreover, it has been mainly developed in an 
Euclidean path integral formulation. In contrast to that, Hamiltonian techniques have re- 
mained less studied. With the Hamiltonian, one can project out the correct ground state by 
evolving an initial wave functional in imaginary time. In continuum theory, some progress 
has been made recently in the non-perturbative regime ^, 0, H]. Accoring to these circum- 
stances there have been only few contacts between lattice QCD and light cone field theory 
(LCFT). 

There is no doubt that LCFT is an important tool for the description of high energy 
interactions. The knowledge of wave functionals in the gauge field configuration space may 
help to calculate light cone wave functions of hadrons. In the following paper we attempt to 
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8, Although the 
1 Of ] offers the advantage 



take advantage of lattice methods in LCFT (for previous work, see 
Hamiltonian is not Lorentz invariant, the light cone Hamiltonian |8j, 
of being boost invariant and has - naively interpreted - a trivial vacuum. On the other hand, 
one would be surprised if QCD looses its non-perturbative vacuum structure in the light cone 
limit. In our opinion much of the complicated vacuum structure of QCD is hidden in the 
constraint equations appearing in light cone QCD. The constraint equations contain zero 
mode solutions which are difficult to solve. These quantum constraint equations have been 
attacked in lower dimensions for scalar theories, but gaug e theories still escape a solution 



111 ] one has been able to solve these 



in higher dimensions. In Nambu Jona Lasinio models 
zero mode equations in the large N c approximation. 

A quantization of scalar light cone field theory on the lattice has been first analysed in 
ref. jsl where also the time coordinate has been discretized. In this reference, special care 
has been devoted to the constraints which arise on the light cone. This approach has not 
found applications. In particular, it is not easily extendable to gauge theories. 

Remarkable progress has been made in light cone QCD with a color dielectric lattice 



theory starting point , ll2j . This approach is based on "fat" links which arise from 
averaging gluon configurations by a block spinning procedure [13I . I4J]. With this method the 
spectrum of glue balls and the pion light cone wave function have been calculated [ijj]. In a 
Lagrangian framework the connection to the original QCD Lagrangian can be easily made, 
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although the numerical accuracy is limited. On the light cone, however, one is prevented 
from approaching the continuum limit, since an effective potential for the link matrices 
M e GL(N) with a non vanishing vacuum expectation value is not allowed. The norm of 
the link matrices M G GL(N), however, should approach unity in the continuum limit. 
This is the reason why we propose to formulate QCD near the light cone. We have 



already analyzed scalar theories [16] and QCD [17], [18| approaching the light cone in a 
tilted near light cone reference system containing a parameter rj ^ parameterizing the 
distance to the light cone. Our work in this paper will follow this idea deriving a lattice 
Hamiltonian which describes the pure gauge sector of QCD and which is suitable for a 
numerical treatment. For QCD, we have already followed the path of maximal gauge fixing 
0, Q] outlined by the Erlangen group in previous works. This way to eliminate 
all gauge degrees of freedom looks very attractive analytically, but numerically it is not 
advantageous. It includes solutions of constraint equations which complicate the form of 
the Hamiltonian. Hence we do not fix the gauge in the following work and try to establish 
a form of the Hamiltonian describing the near light cone dynamics similar to the QCD 
Hamiltonian in an equal time approach, i.e. in terms of unitary matrices describing the gauge 
degrees of freedom and their canonically conjugate momenta. In our lattice prescription, we 
leave near light cone time continuous. It plays a similar role as ordinary Minkowski time, 
therefore, we can follow the conventional method of the transfer matrix in order to derive the 
lattice Hamiltonian from the lattice action. The transversal field strengths are increased in 
magnitude due to the boost into the vicinity of the light cone whereas the longitudinal fields 
remain unchanged. Constraint equations arise in the light cone Hamiltonian framework, 
since the Lagrangian contains the velocities in linear form. The momenta related to these 
velocities obey constraint equations. The constraint equations appear in the near light 
cone Hamiltonian as terms proportional to l/rj 2 . These terms enforce the "equality" of the 
transverse chromo electric and chromo magnetic fields E% = F® k . While the longitudinal 
chromo electric field and the longitudinal chromo magnetic field appear in their usual form 
in the light cone Hamiltonian, the Hamiltonian contains the transverse chromo magnetic 
field squared in an unusual quadratic Z(2) invariant form. The Z(2) invariance, however, is 
broken because the chromo magnetic fields also appear linearly together with the transverse 
chromo electric field. 

The lattice Hamiltonian density depends on an effective constant which represents the 
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product of the anisotropy parameter £ = a_/aj_ and the near light cone parameter 77. If 
one chooses rj — 1 and lets £ — > one obtains a deformed system which is squeezed in the 
spatial (— )-direction, if one uses £ = 1 and lets 77 — > one obtains the light cone limit. 



This equivalence has been advocated before by Verlinde and Verlinde [20| and Arefeva 21]. 
These authors have proposed to implement the strong interaction with such asymmetric 
lattices in order to study high energy scattering, motivating us to proceed in this way. As 
it stands, the (anisotropic) lattice Hamiltonian itself is not usable for Monte Carlo methods 
evolving an arbitrary initial state in imaginary time to the ground state, since the chromo 
electric field strengths i.e. the momenta canonically conjugate to the links appear linearly. 
Therefore we propose to use the translational invariance of the vacuum to add a term l/rfP_ 
in order to cancel the unwanted terms. Naively this amounts to returning to an effective 
lattice Hamiltonian which is proportional to the energy in ordinary Minkowski coordinates. 
For the ground state of the vacuum this seems a reasonable procedure. Applications of the 
light cone coordinates in finite temperature field theory have followed the same route 22]. 
The new effective Hamiltonian contains two parts: The first describing the dynamics of the 
longitudinal chromo electric and chromo magnetic fields is not influenced by the smallness of 
the near light cone parameter rj. The second part containing the transverse chromo electric 
and chromo magnetic terms is enhanced with 77 in the light cone limit. 

We analytically investigate this effective Hamiltonian in the strong and weak coupling 
limit. Such a procedure can direct the search for an appropriate (approximate) guidance 
wave functional needed to improve the convergence of the Hamiltonian Monte Carlo method. 
The strong coupling limit suggests a simple sum of plaquette terms with different weights 
in the purely transverse and (—) -transverse planes. The magnitude of the couplings follows 
the asymmetries existing in the Hamiltonian. In the limit 77 — > the plaquette terms in the 
purely transverse planes are weighted very weakly, i.e. the longitudinal magnetic fields can 
vary freely. The weak coupling approximation identifies the "abelian" fluctuations with their 
modified dispersion relations following the built in anisotropy. It is particularly interesting 
that the light cone limit 77 — > produces long-range correlations in the minus direction, 
which deviate from the local strong coupling ansatz. In fact this anisotropy may help to 
make an ansatz for the ground state wave functional which is especially appropriate in the 
light cone limit. It correlates fluctuations of the longitudinal chromo magnetic fields with 
long strings along the (— )-direction. This ansatz may also point the way to find a solution 



4 



of the quantum constraint of the initial Hamiltonian. 

The outline of the paper is as follows: In Sec. [Til we introduce near light cone coordinates. 
For the sake of clarity, we first establish the methodology in the continuum formulation. 
We derive the continuum Hamiltonian and momentum operator. Furthermore, we motivate 
an effective Hamiltonian making an ordinary Quantum Diffusion Monte Carlo algorithm 
possible. In Sec. II III we switch to the lattice formulation and derive the near light cone 
Hamiltonian from the latticized action with the transfer matrix method. In Sec. IIVI we set 
up the effective Hamiltonian. The time independent Schrodinger equation for the effective 
Hamiltonian is analytically solved for the ground state in the strong and weak coupling limit 
in Sec. El In Sec. |Vl]we variationally optimize an ansatz for the ground state wave functional 
motivated by the strong and weak coupling analysis. It allows to interpolate between these 
two extreme limits and to investigate the rj behavior in the whole coupling range. Finally, 
in Sec. IVIH we present our conclusions and an outlook to future work. 



II. THE CONTINUUM QCD HAMILTONIAN AND MOMENTUM NEAR THE 
LIGHT CONE 

Before we start with the actual derivation of the QCD Hamiltonian and momentum 
near the light cone, we would like to introduce near light cone coordinates similar to the 



coordinates first proposed by |23l . |24j. The transition to near light cone (nlc) coordinates 
might be considered as a two-step process. In the first step, one starts in ordinary Minkowski 
space in the laboratory frame with unprimed coordinates x^ and transforms into a reference 
frame described by primed coordinates x' M which moves with relative velocity (3 along the 
longitudinal direction relative to the laboratory frame. The relative velocity j3 is chosen to 
be given by 



1 H- 77 



The associated Lorentz transformation expressing the primed coordinates in terms of labo- 
ratory frame coordinates reads 

= 1 ( x °-(3x 3 ) 

= 7 (x 3 -(3x°) , 7= (1-/5 2 )" 1/2 . (2) 
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Here x° and x 3 denote the temporal coordinate and the longitudinal spatial coordinate 
respectively in usual Minkowski coordinates. From the boosted frame, one performs an ad- 
ditional linear transformation not included in the Lorentz group which rotates the temporal 
and longitudinal coordinates. It is given by 

1 



x 



X 



(l+rf) x'° + (l-r] 2 ) X 



/3 



x'° - x' 3 



(3) 



Here, x + is defined to be the new time coordinate along which the system evolves and x~ is 
defined to be the new spatial longitudinal coordinate. The transversal coordinates x 1 and 
x 2 remain unchanged. By quantizing a theory on a hypersurface of constant x + , one can 
smoothly interpolate between an equal time quantization and light cone quantization by 
varying the external near light cone parameter 77 from 1 to 0. In the equal time limit 77 = 1, 
the temporal coordinate x + is given by the ordinary Minkowski time coordinate x + = x ' 
and (3 = 0, i.e. the new reference frame is not moving relative to the laboratory frame. In 
the light cone limit 77 — > 0, x + is proportional to the usual temporal light cone coordinate 
and (3 approaches (3 = 1. The nlc energy p + and longitudinal momentum p_ expressed in 
terms of the laboratory energy E and longitudinal momentum p 3 are given by 

1 



P- 



V 



(E-p 3 ) 



VP • (4) 

The second relation in Eq. (j4j) shows that the magnitude of longitudinal momenta is reduced 
by transforming to nlc coordinates. In other words, large longitudinal momenta in the lab 
frame p 3 oc 1/ (a^rj) become accessible by a nlc lattice with longitudinal lattice spacing a_ 
for 77 — ■> 0. This makes nlc coordinates physically very attractive. 

The definition of nlc coordinates Eq. (j3J) induces the following metric: 

/ 1 \ /t7 2 0i\ 
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0-1 
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0-100 
0-10 

\ 1 00/ 



(5) 



with \i, v = +, 1, 2, — , det g = 1. This defines the scalar product 



x y + x y 
x„y + + x + y. 



rfx y 



77 x + y + - x±y± 



(6) 



Note, that the metric has off-diagonal terms which implies that there are terms mixing 
temporal and longitudinal spatial coordinates in the scalar product. This has severe conse- 
quences for a standard Euclidean lattice approach. 

If we put a pair of color charges propagating along the longitudinal coordinate x~ de- 
scribed by a longitudinally extended Wegner- Wilson loop and a stationary target modeled 
by a transversal plaquette at fixed x~ in this reference frame, we can simulate color dipoles 



colliding with a hadron 25j in the light cone limit. In the described way one might be able 



to calculate cross sections between hadrons. For 77 — > 0, we approach the light cone from 



space like distances which is different from the approach of Balitsky 26( who approaches 
the light cone from time like distances closer to scattering experiments. 

For QCD, the pure gluonic part of the Lagrange density in manifestly covariant notation 
is given by 



C = —F^sT^F^ , (7) 



with the non-abelian field strength tensor 

F; v ee d,A a v - d v Al + gr hc A\Al . (8) 

In the following, we restrict ourselves to the color gauge group SU(2) for which the struc- 
ture constants f abc are given by the three-dimensional totally antisymmetric Levi-Cevita 
symbol e abc . By using the nlc metric Eq. §5§ we obtain for the Lagrange density Eq. (J7J) 

2 



k=i ^ ' 



(9) 



Note, that there is a term in the Lagrange density which is only linear in one of the temporal 
field strengths, namely F^ k F^ k . Therefore, the numerical standard approach for lattice 
gauge theory, the Monte Carlo sampling of the Euclidean path integral does not apply for 
nlc coordinates. The reasoning is as follows. If one performs in analogy to equal time theories 
an analytical continuation to imaginary nlc time x + — > — ix^, each temporal field strength is 
replaced by its Euclidean counterpart times an additional factor i. Therefore, the linear term 
yields a complex valued Euclidean action and the integrand of the Euclidean path integral 
is no longer interpretable as a probability density. A similar problem arises for lattice gauge 
theory at finite baryonic densities which is usually referred to as the sign problem. So far, no 
convenient solution has been found. In order to avoid these problems, we stay in Minkowski 



time for the rest of the paper and we switch to a Hamiltonian formulation. We perform a 
Legendre transformation to switch to a Hamiltonian formulation , i.e. we have to express 
the temporal derivatives of the fields by their canonical conjugate momenta in particular 
which are given by the functional derivatives of the Lagrange density with respect to the 
temporal derivative of the correspondent fields: 



K = , 5C \ • (io) 



Therefore, the canonical momenta conjugate to the gauge fields are given by 

k ~ Sd + A' k ~ 6F" +k ~ ~" + ' * k • 

Here, we have chosen the axial gauge A°l = which is quite natural because the temporal 
gauge field A ^ is not dynamical, i.e. there is no temporal derivative appearing in the 
Lagrange function. It acts like a Lagrange multiplier which multiplies the Gauss law G = 
with G given by: 

G = D a _ c Ut + D a k c Ii c k . (12) 

Here D^ c denotes the ordinary covariant derivative - in the adjoint representation - in spatial 
direction \x 

D ; c = d^ + gf^. (13) 

In order to recover the full Lagrangian dynamics, we have to supplement the equations of 
motion by Gauss' law. Hence, the Gauss law has to be imposed as a constraint equation 
on physical states. We express the temporal derivatives of the gauge fields in terms of the 
canonical conjugate momenta by using Eq. (llip . which yields 

d+A a _ = Fl_ = U a _ . (14) 

We may obtain the QCD Hamiltonian and the momentum operator via the energy momen- 
tum tensor, where we have to substitute the temporal derivatives of the gauge fields by the 
corresponding expressions involving the canonical conjugate momenta Eq. (TT41) . If the La- 
grange density for an arbitrary field theory with fields $ r defined by the Lagrangian density 
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£ is a function of the fields itself and derivatives of the fields only, namely C = £($ r , <9 M $ r ), 
the energy momentum tensor in its most general form is given by 

It defines the Hamiltonian density 7i and the longitudinal momentum density V- by 

n = t\, 

V- = T+_ . (16) 
Therefore, for the nlc QCD Lagrangian Eq. ([9]) we find the Hamiltonian density 

(17) 



n = \ n a _rr_ + f»f& + ^ i (u a k - f\)' 

a 



and the longitudinal momentum density 

2 

V- = U a _d^A a _ +J2 u t d - A l ■ ( 18 ) 
fc=i 

This form of the local integrand for the generator V- of longitudinal translations is not 
manifestly gauge invariant. However, if one uses Gauss' law and the definition of the field 
strength tensor one can rewrite V- as 

V- = U a k F^ k + 8 k (U a k A a _) + <9_ (ITM° ) . (19) 

So, the longitudinal momentum density V- may be expressed as a manifestly gauge invari- 
ant object plus some total derivatives along the spatial directions which disappear after 
integration with periodic boundary conditions. We use the symmetrized form 

V- = \ (U a k Fl k + F1AX) . (20) 

The integrated Hamiltonian density H is the generator of nlc "time" translations and the 
integrated longitudinal momentum operator P_ is the generator of spatial translations in 
longitudinal direction: 



H — J d x±dx TC , 

P = I d 2 x ± dx~V- . (21) 



We quantize the theory by choosing the following commutation relations at equal light cone 
time x + : 

[U a m (x),A b n (y)] = -i5 ab 5 mn 5®(x-y), 
[U a m (x),U h n (y)] = 0, 

[A a m (x),A b n (y)} = 0. (22) 

These commutator relations respect the Heisenberg equations of motion. Note, analogously 
for Quantum Mechanics, one has to supplement the Heisenberg equations of motion by the 
Gauss law constraint. In quantum mechanics, the Gauss law constraint translates into a 
restriction of the Hilbert space to the subspace of physical states i.e. states \I/ satisfying the 
Gauss law 

^ 6 lTi(f) +Y^DfU b k (x) S j = Vf.fl. (23) 

Since the Gauss law operator is the generator of gauge transformations the physical subspace 
is given by that part of the entire Hilbert space which is spanned by gauge invariant states. 
The 1 / r/ 2 -term in the Hamiltonian Eq. (ITT)) favors expectation values of transverse chromo 
electric fields and transverse chromo magnetic fields F^ k to be equal in order to have 
a minimal energy. On the other hand, this term introduces terms linear in momentum 
which are difficult to handle, for example with a numerical Quantum Diffusion Monte Carlo 



algorithm [27|, [28|, [29|] which exploits the fact that the time evolution operator is a projector 
onto the ground state when analytically continued to imaginary times. The terms linear in 
the momentum enforce the wave functional to be complex valued in general which spoils 
the whole procedure. These are exactly the same terms which make the nlc action complex 
valued after the Wick rotation in the action based formulation. Hence, the problem reappears 
in the Hamiltonian formulation. However, for Hamiltonian nlc QCD it is possible to define 
an effective Hamiltonian converging to the exact ground state which avoids the problematic 
terms. 

Obviously, the Hamilton operator H in Eq. (1211) is translation invariant and gauge invariant. 
Hence, it commutes with the longitudinal momentum operator P_ and with the Gauss 
operator G: 

[H,P-] = 0, 

[H, G] = . (24) 
10 



Therefore, common eigenstates exist which diagonalize the Hamiltonian and the longitudinal 
momentum operator simultaneously and in addition fulfill the Gauss law. In particular, 
momentum is a good quantum number which is left invariant by time evolution. In order to 
solve the Hamiltonian we are interested in translation-invariant states which are eigenstates 
of the longitudinal momentum operator, i.e. with eigenvalue equal zero. In vacuum, with 
light cone momentum P_ = 0, we can add (I/77 2 ) P- to define an effective Hamiltonian 
density Tl e ff which is only quadratic in momenta: 



H 



eff 



1 £ na _ na _ + F „ ^ + £ 1 ^ + pa _ kp ^ 



k=l 



(25) 



This effective Hamiltonian density is still symmetric under the exchange 



(26) 



but it does not enforce the equality between transverse chromo electric and transverse chromo 
magnetic fields commonly used in the light cone limit also for the quantum field theoretic 
system. 

One finds the ground state |^o) of H by evolving a translation invariant trial state |<&) not 
orthogonal to |^o) with the effective time evolution operator 



|*o) = lim exp [- (H eff - E eff ) r] |$) 



(27) 



related to the effective Hamiltonian. For the details of an explicit implementation of the 
ground state projection operator with a guided Quantum Diffusion Monte Carlo algorithm, 



we refer the reader to 



27 



28 



29]. In order to direct the Monte Carlo into regions of the 



configuration space which have large acceptance rates, i.e. which have a large exact ground 
state probability density for the given configuration, one introduces a guidance wave func- 
tional which is an approximation of the exact ground state. Instead of evolving the trial state 
itself, one evolves a probability density in imaginary time which converges to the product 
of the exact ground state wave functional and the guidance wave functional for asymptotic 
times. Obviously, the application of a guidance wave functional introduces some bias in the 
computation of expectation values. However, in principle one can get rid of this bias by 



applying forward walking techniques 



30 



Hi- 



ll 



The algorithm preserves Gauss' law as long as the guidance/trial wave functional is a func- 
tional of closed loops only. In principle, there are also multiple connected loops possible, 
but then the chromo electric flux must be conserved at each site. In this paper our pri- 
mary objective is to translate the discussed methodology onto the lattice and to determine 
variationally a good starting and guidance wave functional |$) for the Quantum Diffusion 
Monte Carlo evolution based on Eq. (j27j) which is motivated by analytical computations in 
the strong and weak coupling limit. 



III. NEAR LIGHT CONE HAMILTONIAN H ON THE LATTICE 

In order to regularize the continuum formulation we go over to the lattice formulation. 

n 

In a previous paper [32J we have started deriving the Hamiltonian for gauge theories on the 
lattice from the path integral via the transition a + — > 0. Here, we go through the procedure 
in detail. We introduce in four- dimensional space SU(N) link variables Ui(x) connecting 
the site x with the 4D site x + ei (e$, £ = +,—, 1, 2, is a unit vector in 4D space-time) in the 
following way: 

(x+ei \ 
ig J dy»A«(y)^\ , (28) 

where V is implementing path ordering from left to right with increasing y^ and a a are 
hermitian color generators. In the following, we restrict to SU(2) where the a a are given by 
the Pauli matrices. The hermitian conjugate of the link variables, Uj(x), connect the site 
x + ei with the site x in reverse order. Plaquettes Uij(x) are related to the field strengths 
Fij(x) and have the usual form 

U tJ (x) = Ui^Ujix + e t )Uj(x + ej)U](x) . (29) 

Expanding a plaquette around its center x + % + % in orders of the lattice spacing one 
obtains 

U tJ (x) = I + iga iaj F^ - ^a^F^^ + . . . . (30) 

Here a, denotes the lattice spacing for direction i, i.e. we allow in general for different 
lattice spacings in the temporal, longitudinal and transversal directions. A correspondent 
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expansion is obtained for UjAx). Therefore, the sum over color indices a = b of a product 



of field strengths is given in the limit cij 

4 



-> as follows 
Tr 



1 -Ref Uy{x) 



Or more general 



-Tr 



g'aidjdkai 
Here, Re (U) and Im (U) are defined as 

U + W 



Im( Uij(x) 



Re (U) 



lm(U) 



Im( Ukiix 



U-W 



(31) 



(32) 



(33) 



2 ' v ' 2i 

In Eq. ( 152"]) the two plaquettes t/y and begin and end at the common site x. Note, that 
Eq. (131]) and Eq. (132]) are representations of the field strengths squared terms which are 
only valid in leading order of the lattice spacing. So far, there is no improvement included. 
By using the relations Eq. (13~T]) and Eq. ( 132]) . we may rewrite the continuum nlc Lagrange 
density Eq. ([9]) in terms of plaquettes such that it is recovered in the naive continuum limit. 



Similar to the equal time case 



33j, one can fix inside the path integral on the lattice a 



maximal tree of links to arbitrary group elements. A maximal tree of links is a tree to which 
no more links can be added without forming a loop. By doing so, the path integral itself and 
expectation values of gauge invariant operators are not affected. Hence, we fix all time-like 
links U + (x) to U + (x) = 1 in the following. This corresponds to the temporal gauge A + = 
and one obtains for the lattice analog Si at of the action S = J d A x C: 



lat 



X 

a + a. 



a_i_a_ 



Tr 



1 - Re U-(x + e+)W_(x) 



Tr 



1 - Re Unix 



+ 



Tr [lm(u k (x + e + )Ui(x)^) Im(V_ fc ( , ) 



+ - 



v 2^ 

k 



Tr 



1 - Re U k (x + e+)Ul(x) 




(34) 



Therefore, the QCD path integral on the lattice in the 
Haar measure dU is given by 



n n du > 

L x 3=1,2,- 



jSi, 



gauge and with the SU(2) 



(35) 
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In order to obtain the lattice Hamiltonian, we would like to go over from the action-based 
path-integral to a Hilbert-space formulation of the near light cone QCD lattice gauge theory 
in the following, letting the time-like lattice constant approach zero. The method is similar 
to the transition from the action to the Hamiltonian in ordinary Euclidean SU(2) lattice 
gauge theory carried out by Creutz [33 1. 

The procedure consists of two steps. First, we construct the transfer matrix T. Second, we 
define the space on which it acts and rewrite the transfer matrix in terms of the conjugated 
momenta of the links and extract the lattice Hamiltonian by identifying the transfer matrix 
with the time evolution operator which propagates the system from one time slice to the 
next. 

Note, that the lattice action Eq. (|35p is local in the temporal direction. Each piece is 
connecting two adjacent time slices x' + = x + + a + and x + which means that the path integral 



factorizes into a product of transfer matrices T(x 

2 6 



,+ ,x+) 



T 



II exp 



1_ 2 

g z a + a 



Tr 



nf 2 2 a_ 
exp <^ i — 7] — Tr 

. x,k 



exp < i — Tr 

T 



J^exp 



— i 



1-Re( [/_(£, x' + )W_(x,x 
I - Re(u k (x, x' + )Ul(x, x + ) ) 
Im( Uk(x, x' + )Ul(x, x + ) \ Im( U-k(x, x + 
1 - Ref U 12 {x,x + 



x 



2 a + a 



x 



g 2 a\ 



Tr 



(36) 



Here, x denotes a lattice vector in the three dimensional spatial sub lattice. If we denote by 
the set of links U(x + ) an entire spatial lattice configuration at time x + , the transfer matrix 
T evolves the configuration U(x + ) at the time slice x + to the configuration U(x' + ) at the 
neighboring time slice x' + in our convention. The construction of the Hilbert space and the 
transcription of the temporal plaquettes in terms of momenta canonically conjugate to the 
links is similar to the steps performed in [33]. The interested reader may find the explicit 



14 



calculation in appendix [B] One finally obtains for the lattice Hamiltonian 



Hi 



lat 



2 a_ 



^-E A 

k,a 1 



IIJ(x) - - Tr 




g 2 a_ v^™ /_*•> 2 a_ 



„ 9 2 cl\ 



Tr 



^ Imf £/_ 
2 V 

1 -Re( U 12 (x 



(37) 



Here, the operators are canonically conjugate to the link operators and they obey the 

following commutation relations 



[nj(a),fy(f)] 
n|(f),c/j,(f)" 
[n-(f),n5,(f)] 
[n^f) 2 ,!!^ 



0", 



-U}{x) y , 

i e afec n^f) ^ , 
o . 



(38) 



In analogy to the continuum Hamiltonian density cf. Eq. ffT7|) we introduce the lattice 
Hamiltonian density 

H\a,t 1 Hat 



lat 



(39) 



Mat ALiV 2 

The lattice anisotropy parameter £ is given by the ratio of the longitudinal lattice spacing 
to the transversal lattice spacing 



a_L 



(40) 



Furthermore, in order to simplify the notation, we have introduced the coupling constant A 
which is related to the ordinary SU{2) lattice gauge theory coupling (3 by 

2 

R = - 



g 2 



(41) 



Therefore, we obtain for the Hamiltonian density on the lattice 



U 



lat 



1 ' 2 E^E + ^ ATr 



y 1 - — 



1 - Re Unix 



Il a k (x) - V\ Tr y y Im( C/_ fc (f) 



(42) 
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One observes that the energy density H\ at only depends on the effective constant fj defined 
as the product of the anisotropy parameter £ = a_/a_i_ with 77 instead of both of them 
separately 

V = £ • V • (43) 

Very clearly one can vary two independent parameters X = 4/g 4 and fj. The fj variation may 
be interpreted in two parametrically distinct but physically equivalent ways. If one chooses 
77 = 1 and varies £, one simulates an effective equal time theory with a ratio of lattice con- 
stants £ = a_/aj_. In the limit £ — > one ends up with a system, which is contracted in the 



longitudinal direction. Verlinde and Verlinde 20( and Arefeva 2l| have advocated such a 
set-up to describe high energy scattering. A contracted longitudinal system means that the 
minimal momenta become high in longitudinal direction and this looks a promising starting 
point for high energy scattering. It is obvious that this limit leads to the same physics as the 
limit 77 — > and £ = 1, i.e. as the light cone limit with equal lattice constants in longitudinal 
and transverse directions. 

In both limiting cases, i.e. for fj — > 0, the near light cone Hamiltonian is dominated by 
the term proportional to (1/fj 2 ). Therefore, in the light cone limit the transverse chromo 
electric fields 11^ should become equal to the scaled transverse chromo magnetic fields 
Tr[cr a /2 Im([/_fc)]. This is a form of electro-magnetic duality characteristic of light cone 
gauge field theory. 

In the following we set £ = 1 bearing in mind that the physical ratio of longitudinal to 
transverse lattice spacings for 77 ^ may be modified by quantum corrections from the 
QCD dynamics. 



IV. EFFECTIVE NEAR LIGHT CONE LATTICE HAMILTONIAN 

To obtain the same cancellation of linear terms in 11^ in the effective lattice Hamiltonian 
as in the continuum Eq. ( 1251) in order to make a guided Diffusion Monte Carlo in principle 
possible, we add P-,i a t to the lattice Hamiltonian density 

W e ff,lat = T^lat H ^-.lat • (44) 
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Here the density P_i a t is defined as 

P-,lat 



V. 



,lat 



Mat 
1 



1 1 



E( n ^ f )" Tr [y Im (^( f )) 



x,k,a 



+ Tr 



- Im^f/„ fc (f) 



(45) 



In the naive continuum limit, i.e. for infinitesimal a_ Eq. (j4"5T) becomes the generator of 
translations along the longitudinal direction. However, P_ iat does not generate translations 
on the lattice for finite lattice spacings. As a consequence, translation invariant states on 
the lattice are not exact eigenstates of P-jat- There are higher order corrections in a_ which 
prevent P_ ; i a t from being the exact longitudinal lattice translation operator. 

In a numerical simulation with an explicit implementation of the ground state projection 
operator one has to ensure that the substitution of the lattice Hamiltonian by the effective 
lattice Hamiltonian is justified. If the relative magnitude of the corrections with respect 
to the ground state energy is of the order of the time evolution step size At of the Quan- 
tum Diffusion Monte Carlo, then the induced defect in the time evolution is effectively of 
quadratic order in At. Hence it can be safely neglected due to the fact that the Quantum 
Diffusion Monte Carlo algorithm itself is only valid up to quadratic order in At. In order to 
quantify the quality of the substitution it is important to measure the typical magnitude of 
the fluctuations (Pf lat ) of P-^t around its expectation value (P_ 5 i at ) where the expectation 
values are computed with respect to the translation invariant trial/guidance wave functional 
to which the projection operator is applied. Both expectation values are equal to zero for 
the exact generator of longitudinal translations. This is not true for P_ 5 i at . In order to 
minimize the defect in the time evolution, the expectation value of P_,i a t with respect to the 
trial wave functional has to be equal to zero and its relative fluctuations have to be of the 
order of the time evolution step size as discussed. Therefore, the trial wave functional has 
to be selected accordingly. 

The effective lattice Hamiltonian can then be chosen as 



n 



eff.lat 



k.a 



1 2 £ U £ n " 



x) 



+ - A Tr 



1 - Re Un(x 



U a k (x) 2 + A 



Tr 



On 



— Im( U-k(x) 



(46) 
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By construction, the effective lattice Hamiltonian Eq. (|4"6]) is equivalent to a naively latticized 
version of the effective continuum Hamiltonian Eq. (1251) . For fj = 1 this effective lattice 
Hamiltonian is very similar to the traditional Hamiltonian used in equal time lattice theory. 
They differ in the potential energy terms for the plaquettes. Instead of the usual 
Tr[l — Re(E/_fc)] term resembling the field strength squared in the naive continuum limit, 
the effective nlc Hamiltonian has the form (Tr[cr a /2 Im([/_fc)]) 2 which corresponds to the 
plaquette in the adjoint representation. These terms which coincide in the continuum limit 
have different finite lattice spacing corrections. 

Note that the effective Hamiltonian Eq. (HSl) has a symmetry which the original Hamiltonian 
Eq. (1421) did not have, namely it is invariant under a Z(2) transformation of the following 
form 

Uk(x±, x~) — > z Uk(x±, x~) V x± and x~ fixed , z G Z{2) . (47) 
Under this transformation, the longitudinal-transversal plaquettes 

U-k(x±, x') and U^k(x±,x~ — 1) involving transversal links belonging to the longitudinal 
slice x~ transform like 

U-k(x±,X~) -> 2(/_ k (fi,^) 

U-k(x±,x~ - 1) — >■ z U-k(x±,x~ - 1) . (48) 

Of course, this symmetry can be spontaneously broken. In order to preserve the symmetry 
properties of the original Hamiltonian we have to restrict ourselves to the phase in which 
the symmetry is spontaneously broken. The order parameter of the phase transition is the 
expectation value of Tr Re In the symmetric phase, the expectation value is equal to 
zero and in the broken phase it acquires a non-vanishing expectation value 



Tr 



Re U. 



= Z(2) symmetric phase 
7^ Z(2) broken phase 



The light cone limit fj — > enhances the importance of transverse chromo electric and 
magnetic fields similar to the full nlc Hamiltonian without the unwanted linear terms in the 
momenta. The resulting vacuum solution should be a plausible extrapolation of the vacuum 
solution of QCD. 
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V. ANALYTICAL SOLUTIONS OF THE EFFECTIVE LATTICE HAMILTO- 
NIAN 

With regard to a subsequent implementation of a guided diffusion Monte Carlo it is 
important to know as much as possible about the true ground state. Analytical solutions of 
the effective lattice Hamiltonian are possible in certain regions of the parameter space given 
by (A, 77). In particular, we would like to analyze the behavior of the ground state wave 
functional, i.e. the vacuum state, when the effective parameter 77 — > makes the vacuum 
approach the light cone vacuum. Therefore, we have a closer look on the strong (A << 1) 
and weak coupling (A >> 1) solution of the Schrodinger equation for the effective lattice 
Hamiltonian in the following. 



A. The strong coupling solution of the effective lattice Hamiltonian 

In this section we investigate the strong coupling limit of the Schrodinger equation for 
which we are able to find analytic solutions. In the strong coupling limit g » 1, i.e. A << 1 
the effective Hamiltonian density Eq. (I46p is dominated by chromo electric fields which 
represent the kinetic energy terms. In comparison with the kinetic energy, the potential 
energy terms are suppressed by factors of A = 4/g 4 . Therefore, we may interpret the 
effective Hamiltonian density as an unperturbed part T plus a small perturbation A V po t 



lat 



T+AV, 



pot- 



(50) 



Here the kinetic energy density T is given by 



T 



1 2 



ALiV| 4 v/A 



E 



(51) 



The potential energy density A V pot represents a small perturbation 



pot 



I 1 2 ^111 



N.Nl 4 v/A 



1-l-Tr 



Re UJx) 



+ 



1 - -Tr 

2 



Re Un(x 



(52) 
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In order to write the potential energy density A V po t in the given form Eq. (1521) . we have 
used the following identity 

2 /i \ 2 



<7 



Im U-k(x) 



1 



-Tr 



Re I7_ fc (£) 



(53) 



We perform perturbation theory in A. Then the ground state |\l/o) as wen as the ground 
state energy density e are written as a power series in A where "(n)" denotes the n-th order 
correction 



i* > = E An l* 



(n) 




e 



n=0 

oo 

E A ' 

n=0 



» 



(54) 



The unperturbed Hamiltonian T is a sum of quantum rigid rotators, one for each lattice site 
and for each spatial direction 34]. The spectrum of each Yl a ^ a2 * s given by E\ = 1(1 + 1) 



with I G (0, 1/2, 1, . . .) in SU(2). Each eigenvalue E\ is (2/ + l) 2 -fold degenerate. Therefore, 
the unperturbed ground state l^ ^) of T is the state which has I = for each rotator. It is 
annihilated by all the momentum operators 



W 3 (x) 



(o) 



W x,a A V j G {1, 2, — } . 



(55) 



This state does not depend on the Uj(x) in the link-coordinate representation, i.e. is a 
constant and is non-degenerate. The corresponding ground state energy is given by 



^ = o. 



(56) 



The space of states may be constructed from the ground state I^q ^) by applying the link 
operator in a given representation (/) which is then again an eigenstate of Y^II|(x) 2 with 

a 

eigenvalue Ei 



(57) 



Note that the representation index (I) of the link explicitly refers to its S'L r (2)-representation 
whereas links without a representation index are defined to be in the fundamental represen- 
tation 



U,{x) = Uf /2 \x) 
20 



(58) 



Due to the non-degenerate ground state we may apply standard Raleigh- Schrodinger per- 
turbation theory. In general, the first order correction to the ground state reads 

1 



(i) 



c(°) _ T 



-V, 



pot 



(0) 



(0) 



(59) 



The correspondent first order correction to the ground state energy density is given by 



(0) 



V, 



pot 



(0) 



(60) 



It is a Haar integral over the whole configuration space which is given by 

xx) - 1 V pot {U)\{dU, 



X.J 

112/3 



rf a\ 



This yields a total ground state energy density in the strong coupling limit 



1 1 



+ 2fj 2 )VX + 0(X 3 / 2 ) 



(61) 



(62) 



In order to compute Eq. ( 1591) we use the fact that the trace of the plaquette U^ix) and 
the squared trace of the plaquette LLfc(x) minus one are eigenstates of the kinetic energy 
operator with eigenvalues t_ and t±, respectively (cf. Eq. (1A2|) and Eq. (1A3I) in appendix 

EI) 



T Tr 



Re U»(x 



T 



Tr 



Re tf_ fc (x) 



t_ Tr 



Re Un(x 



(o) 



o 



Tr 



Re U^kix) 



(0) 



(63) 



The eigenvalues t_ and t± are given by 

1 12 2 



_N_Nla\yf\f) 2 \ 4 

112/ 1 
1 + - 



2 . 



(64) 



_iV_iV| 4 ^ V V 2 , 
The factor 3/4 in t- is related to the fundamental representation (/ = 1/2) of the plaquette 
and the factor of 2 in t± arises from the squared trace of the plaquette minus one in the 
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fundamental representation which is equivalent to the trace of the plaquette in the adjoint 
representation (1 = 1). Hence, the first order correction to the ground state wave functional 
is given by 



(i) 



^j^rf Tr[Re(f/ 12 (£ 

X 

5>' 



16 1+^) 



Re Udx) 



(0) 



(65) 



The state \^q) does not contain any products of plaquettes involving field strengths at 
different spatial positions. Therefore, to this order in perturbation theory, the ground state 
wave functional factorizes in a product of single plaquette wave functionals similar to the 
vacuum wave functional obtained for an equal time lattice Hamiltonian [281] 



l*o> 



1 



1 1 



16 1 



rf 



-f? 2 Tr 
3 ' 



J2 ( Tr 



Re U 12 (x) 



Re U- h (x) 



(0) 



H ex Pi o AfTr [Re([/ 12 (£) 



(66) 



A 



+ ^V(Tr RelUJx) 

16 1 + fj 2 ^ V V * v ' 



(0) 



0(X 2 



In the wave functional, the purely transversal plaquettes U12 involving the longitudinal 
chromo magnetic fields are suppressed by rj 2 in the light cone limit fj — > 0. To this order in 
perturbation theory, the strong coupling ground state wave functional Eq. fl66|) respects the 
Z(2) symmetry of the effective Hamiltonian which the full Hamiltonian, however, does not 
share. 



B. Weak coupling solution of the effective lattice Hamiltonian 

In the weak coupling regime, i.e. g —>■ or A — > 00 the effective lattice Hamiltonian 
Eq. ( H6l) in SU(2) depends on a triplet of free U(l) gauge fields and their corresponding 
momenta. To reduce the Hamiltonian into this form it is convenient to substitute the gauge 
field gA?(x) in Eq. (fJED by a rescaled gauge field A*(x) (cf. Eq. (joTjl ). Note that all 
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vector indices throughout this section refer to a flat space metric equal to the unit matrix. 
Furthermore, is the totally antisymmetric Levi-Cevita symbol with eu- = 1. In the 
g — > limit, the field strength tensor reduces to the chromo magnetic field B£(x), which is 
the i-th spatial component of the lattice curl of A a (x) and which is rescaled to Bf(x) 



ge Um [A a m {x) 
W(x) 



1,2, 



el)] 



gA?(x) 

B a (x) 

= (67) 

Similarly to the equal time theory 27J, |28[] one can expand the effective lattice Hamiltonian 
in a power series in A . The expansion of the potential energy is straightforward. The 
kinetic energy of the effective lattice Hamiltonian is a sum of the Casimir operators acting 
on SU(2). Each of them represents a Laplace-Beltrami operator on the curved manifold of 
SU(2). The expansion in a power series of this operator yields in leading order a flat space 
Laplacian in three dimensions given by 

Vm(f) 2 = -JZ^S— 
a 3 r^(^) 2 



(68) 



Hence, the n a ,v4 a obey effectively the following commutation relations 



n?(aO,n}(0 




0. 



(69) 



The described expansion of the effective lattice Hamiltonian in the weak coupling limit yields 
in leading order 

1 I I 



H, 



effjat 



iV_iV 2 a 4 



-^{aH^x) 2 + \B^{xf 



+ o 



^2 A 5/4 



(70) 



This Hamiltonian is equivalent to the abelian limit and the order A 5//4 corrections represent 
the triple gluon vertex g AAA. 
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Instead of solving the ground state in terms of the gauge variant fields as done in ref. 



351 ]. we express the kinetic energy operator acting on the gauge fields in terms of effective 
operators which act on chromo magnetic fields Bk- These are gauge invariant in the abelian 
limit. By doing so, we obtain a ground state wave functional which depends only on gauge 
invariant objects. This avoids an otherwise necessary projection onto a gauge invariant 
subspace of the Hilbert space. After transforming the Hamiltonian into Fourier space, several 
unitary transformations convert the Hamiltonian into a Hamiltonian of decoupled harmonic 
oscillators. The necessary unitary transformations are similar to transformations performed 



for a compact equal time U(l) Hamiltonian in ref. 35[. However, additional factors due to 
the nlc metric appear which can be traced in the computation. Once the harmonic oscillator 
Hamiltonian is obtained, the ground state wave functional and the ground state energy e 
are easily found 



16 1 

k 



2 2 2 _ ^ 

T] sm ( — j + T] sin I — j +sinly 



(71) 



Here L denote the lattice momentum values 



27T 

hi = —Hi m = 0, . . . , Ni - 1 . (72) 

In Fig. CD we show the dimensionless energy density Eq. ( 17T|) for a 16 3 -lattice as a function 
of rj. A leading 1/rj 2 -dependence of the effective energy density is obvious from Eq. (17T|) 
and arises from the 1/fj dependence of the light cone energy and the fj dependence of the 
volume V = N"j_N_a\a_. This dependence is scaled out in the figure. In the abelian limit, 
the energy density is given by the dispersion relation summed over all modes, times the 
color degeneracy factor. If we identify pi = sin(7r Tii/Ni) with the latticized version of the 
i-th component of the gluon momentum p^, then the nlc dispersion relation u n i c of a free 
gluon gas is given by (cf. Eq. (jlj)) 



^ic = - (\Jp\+pI + pI-pz 



(73) 

P3=P-/V 



Here ps refers to the longitudinal mode in the laboratory frame and p- refers to the longitu- 
dinal mode in the nlc frame. Hence, by summing uj n i c over all modes and taking into account 
that the total longitudinal momentum adds up to P_ = we obtain Eq. (1711) . The ground 
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FIG. 1: Rescaled dimensionless energy density eoa^fj 2 of the effective nlc Hamiltonian in leading 
order of the weak coupling limit for a 16 3 -lattice as a function of fj 

state wave functional is a multivariate Gaussian wave functional in the chromo magnetic 
fields where T % ~ (x — x') denote the matrix elements of the covariance matrix 



x,x' a,i,j 



v& = exp^-^^^|sf(f)rf(f-f')|s«(x') 
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Jfj(x - f) 
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V 



(74) 





7 <j (f-f) 
ff^^x — x') J 

Here 7^ denotes the spatial part of the covariance matrix. It depends only on the relative 
distance x — x' of the chromo magnetic fields in the wave functional 

^ 1 1 



7^2 - X 



' 2 N_N\ 

[f sin {h/2) 2 + fj 2 sin {k 2 /2) 2 + sin (k^/2) 2 ] ~ 1/2 



(75) 



The function 7^(x — x') is real due to the invariance under space reflections. In Fig. [2] we 
show 7*7(0) for a 16 3 -lattice as a function of 77. The asymptotic behavior of 7^(0) in the 
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FIG. 2: Diagonal element 7«(0) of the spatial part of the covariance matrix for a 16 3 -lattice as a 
function of fj (solid line). Its asymptotic behavior in the light cone limit fj — > 0, 7^(0) ~ 0.038/r; is 
shown by the dashed line. 



light cone limit 77 — > can be computed by summing all modes with k 7^ and k. 
Eq. (J75l) . For a 16 3 -lattice, it is given by 

0.038 



7«(0) 



7/ 
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. 



in 



(76) 



For fj — 1 the 3x3 covariance matrix Eq. (17^1) equals the covariance matrix which was 
found by Chin et al. 27J, [28( for an equal time theory since our Hamiltonian coincides with 
the correspondent equal time Hamiltonian in the weak coupling limit. For small values of 
fj the chromo magnetic field in the longitudinal direction B a _ oc F® 2 is suppressed in the 
wave functional by a factor of fj in comparison with the other field strengths. On the other 
hand, the chromo magnetic fields in transversal directions B± oc F_2 and I?f oc F_x are not 
suppressed. 

We compare correlation matrix elements 7^(2: — x') for Ax 7^ with the matrix element 
at Ax = by forming the ratio R(Ax) 

-fij(Ax) 



R(Ax) = 



7*(0) 



(77) 
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In Fig. [3j R(Ax) is shown for a 16xl6-lattice and for three different values of fj, namely 
fj — l,fj — 10 _1 and fj = 10 -2 . For reasons of presentability, we have restricted ourselves 
to a 2-dimensional section through the 3-dimensional lattice spanned by x± = x 1 and x~ 
at x 2 = 0. This representation allows to see the anisotropy developing for very small fj. 
Here and in the following, the notion "off-diagonal in position space" refers to Ax ^ 
whereas "diagonal in position space" refers to Ax = 0. For fj = 1, the covariance matrix 
has only weakly-off diagonal contributions in position space. Therefore, it is reasonable to 
consider the weak coupling wave functional in the diagonal approximation 28) as a product 
of single plaquette functionals. For decreasing fj one observes that the correlations among 
plaquettes separated along the longitudinal direction become more and more important. In 
the light cone limit, every plaquette is equally correlated with any other plaquette which is 
longitudinally separated from the first one. 

However, for not too small values of fj, at least an effective description by a product of 
single plaquette wave functionals is possible. In Sec. IVHI we discuss a possibility to include 
long range correlations in the wave functional by a combined optimization and Quantum 
Diffusion Monte Carlo method. 
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f] = 0.01 



FIG. 3: Ratio of covariance matrix elements R(Ax) as a function of the separation Ax for a 2 
dimensional 16x16 lattice at three different values of r] 2 . 
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VI. VARIATIONAL OPTIMIZATION OF THE GROUND STATE WAVE FUNC- 
TIONAL 



In the last two sections we have analyzed the strong and weak coupling behavior of the 
Hamiltonian and its ground state. We have seen that in the strong coupling limit the ground 
state wave functional may be approximated by a product of single site wave functionals. 
Also in the weak coupling limit for not too small fj the bilocality of the chromo magnetic 
field strength is less important. In the following we construct an effective wave functional 
which smoothly interpolates between the strong and weak coupling solution. In addition we 
would like to choose the ground state wave functional in such a way that it is not invariant 
under the unwanted additional Z(2) symmetry of the effective Hamiltonian in which the 
linear momentum terms are compensated by the translation operator. Therefore, me make 
a variational ansatz of the ground state wave functional for the whole coupling range which 
contains a product of single site plaquettes with two variational parameters p and 5. We 
denote the normalization constant by N 

2 



MP,$) = iVjJexp^pTr [Re(c/_ fc (x) ) ] +5Tr [lle([/ 12 (£)) 



. k=i 



(78) 



With this normalized wave functional we variationally optimize the energy expectation value 
€o(p, 5) of the effective Hamiltonian which is given in terms of plaquette expectation values. 
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(79) 



The explicit dependence of the energy expectation value on p and 5 comes from the kinetic 
energy terms in Ti e ff- There is an implicit dependence in the plaquette expectation values 
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which are computed as averages over link configurations generated by the probability density 

dP{u) = \Mp,s)\ 2 i[ vu ^ ■ ( 8 °) 

With the special form of our trial ground state wave functional Eq. (ITS]) , the energy ex- 
pectation value of the effective Hamiltonian Eq. f|46l) coincides with the energy expectation 
value of the full Hamiltonian Eq. (H2l) . Even if we do not use the invariance of the trial 
wave functional under translations, the expectation value of the longitudinal momentum 
operator with respect to the trial wave functional Eq. (175|) vanishes identically. This is due 
to the fact that the expectation value of the chromo electric field operator II" (y) times an 
arbitrary functional G({U}) of the links with respect to a purely real valued exponential 
wave functional obeys 

(tfol IW) G({U}) |* ) = - (* | G({U}) W 3 {y) |* ) . (81) 

The above relation Eq. (181 p may be interpreted as a "partial" integration rule and is proven 
in appendix |A] Hence, the ground state wave functional Eq. (ITS!) minimizing the energy 
density Eq. (179]) optimizes simultaneously the effective and the full Hamiltonian. In order to 
optimize the ground state wave functional we sample the probability distribution Eq. (180]) 



with a local heat bath algorithm 36j on a 16 3 -lattice and measure the expectation values of 



the plaquettes and the squared plaquettes with the bootstrap method 37] using an initial 
sample size of 500 and a bootstrap sample size of 1000. We compute these expectation values 
as a function of the parameters p and 5 on a 50 x 50 grid where each of the parameters varies 
in the interval [0, 10] with a sterilize of 0.2. This yields a set of 2500 distinct expectation 
values which we interpolate with polynomials of fifth order. For a first coarse estimate of 
the optimized parameters, we find the minimum of Eq. (1791) with a standard Mathematica 
minimization routine. 

For the fine determination of the optimal parameters we then generate 50 different pairs 
with energy expectation values less than three percent higher than the energy at the coarse 
estimate of po, 6o- Finally, we fit these energy values with a quadratic form Eq. ( 1821) centered 
at the optimal values (p , ^o) where the linear term in the taylor series vanishes due to the 
minimum condition 

Cs Ct i I »-»° 




d s d p d s ds ) \ 
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The described method is tested by comparing our results with the variational results of 



Chin et al. [28|] who optimized a one parameter (2 p = 2 5 = achin) wave functional of 
the form given in Eq. (I78p with respect to the standard equal time Hamiltonian containing 
only plaquette terms without anisotropy. Note that Chin's results have been obtained on a 
4 3 -lattice, but the authors show that the dependence of the energy density and the optimal 
wave functional parameter on the lattice size is small. We find 0.5% agreement between the 



28]. 



results of our method and the values obtained by Chin et al. 

Next we apply the described optimization method to the nlc Hamiltonian. The optimized 
energy density is presented in Fig. H]as a function of A for different values of fj 2 . The 1/fj 2 
divergence is scaled out. The curve has a a/A behavior for strong coupling and is independent 
of A for weak coupling as found in Sec. [V] 
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(83) 



In Figs. [5] and [6], we present the variationally optimized wave functional parameters po an d 
5q as a function of A for different values of fj 2 . The parameters are divided by a factor y/\ 
such that they become constant in the weak coupling limit (A — > oo). The uncertainties on 
the variational parameters are typically 5% and are larger in the region where the Hamilto- 
nian with the adjoint plaquette in (— A;)-direction induces a phase transition. Therefore, in 
principle only couplings in the weak coupling region above A = 7 are meaningful where the 
Z(2) symmetry is spontaneously broken. 

By using the strong coupling solution from Eq. fl66l) and the diagonal part of the co- 
variance matrix at Ax = of the weak coupling solution Eq. (|74l) . we get analytically the 
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following estimates for p and 8q 








^7^(0) for A >> 1 

| A fj 2 for A << 1 

V\ f] 2 7^(0) for A >> 1 



0.038/77 for fj -> 
0.454 for fj — ► 1 



for A << 1 



(84) 



The variationally determined parameters are in good agreement with the analytic predictions 
in the strong coupling regime. In the weak coupling regime the optimal parameters differ 
from the analytical estimates Eq. ( 1B4"|) . In both cases the analytic predictions disagree 
more for small fj. This is natural, since the light cone limit fj — > builds up correlations 
among plaquettes separated along the longitudinal direction. The parameters optimizing 
our product of single plaquette wave functionals effectively describe these correlations and 
adopt values which differ from the weak coupling estimate given by the diagonal entries of 
the covariance matrix Eq. (154"j) . 

In the following we analyze the fj dependence of the optimal wave functional parameters 
for fixed values of A which lie in the physical relevant region above A = 7. We show in Figs. [7] 
and [5J the optimal wave functional parameters po,^o divided by \f\ 71 (0), i.e. V'A 7^(0) for 
fj = 1, which is the expected behavior for the equal time Hamiltonian. This way we can 
show the variations of the wave functional parameters in the light cone limit. For a direct 
comparison, we plot the analytical weak coupling prediction Eq. ( I84l) by dotted lines in 
the same figures. The analytical results for p (Eq. (1541) ) overestimate the variationally 
determined values, whereas the analytical predictions for 5 (Eq. (1541) ) underestimate the 
optimized parameters as a function of fj. Here again, the large difference for small fj originates 
from the effective description of long range correlations by the parameters of our ground state 
wave functional in this parameter region. For sufficiently large values of A, the fj behavior for 
Po and £0 becomes universal and independent of A. We determine functions f p and fg which 
describe the deviations of the variationally optimized wave functional parameters from the 
the weak coupling limit \f\ 71 (0) at fj = 1 (cf. Figs. [7]and[E]) 



<5 (A,r?) 



^71(0) f P (X,fj) 

Va 7i (3)/*(a, fj). 



(85) 
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In the extreme weak coupling limit A — > oo and close to fj — > 1, each of the functions f p and 
fs may be described by linear functions of fj. Therefore, it is reasonable to assume that f p 
and fs can be approximated by expansions around A — > oo and fj = 1 



fi(X,Tj) = C 

i = p, 5 . 



A 



c 5,; 



Z\2 



fj) 



(86) 



The coefficients Co^ represent the effective single plaquette equal time wave functional pa- 
rameters. A good fit of the parameters c^j, ■■■,c^ > ^ minimizing y 2 in the range A G [10,95] 
and fj G [0.15, 1] gives the coefficients tabulated in Tabled! This analytical parameterization 



i 


CQ,i 


Cl,i 


C2,i 


C3,i 


C4,i 


C5,i 


P 


0.90 


-1.74 


0.72 


4.06 


-0.40 


-0.14 


5 


0.95 


0.93 


-1.21 


-3.22 


-0.83 


0.32 



TABLE I: Coefficients of Eq. (j86j) obtained from least square minimization. 



of the ground state wave functional allows to smoothly interpolate between ground state 
wave functionals belonging to different coupling constants and different values of fj in the 
physical relevant coupling constant region. Furthermore, the given form induces generically 
a vanishing expectation value of P_ ; / a t which makes it optimal for the use in a guided diffu- 
sion Monte Carlo as discussed in Sec. IIV1 Since it is an approximation to the exact ground 
state it may be used for further qualitative investigations: In a forthcoming paper we plan 
to determine hadronic cross sections by simulating how a color dipole moving along the light 
cone hits a neutral hadron localized at x~ = 0. With the parameterization of Eq. fl8"B"|) we 
are able to extrapolate the parameters of the wave functional to fj = 

0.87 1.65 v 



Po(A,0) 
*o(A,0) 



0.65 - 
0.05 



A + A 2 
0.04 1.39 



A 



A 2 



Vx 
Vx. 



(87) 



At fj — 0, the color dipole can be represented by a longitudinal-transversal Wilson loop 
extended in x~ direction and the simplified target can be modeled by a transverse plaquette. 
Varying the impact parameter one can sample the correlation function of the two gauge 
invariant objects and thereby obtain the profile function. 



33 




A 

FIG. 4: Optimized energy density as a function of A obtained from the simulation on a 16 3 lattice for 
three different values of r/ 2 . The red shaded area corresponds to the phase transition region for all 
values of fj 2 . The dotted lines show the predicted analytical strong coupling behavior. The arrows 
indicate the expected asymptotic behavior for weak coupling which is a constant independent of 
A. 
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FIG. 5: Optimal wave functional parameter po(\,fj) as a function of A obtained from the simu- 
lation on a 16 3 lattice for three different values of fj 1 . The red shaded area corresponds to the 
phase transition region for all values of fj 1 . The dotted lines show the predicted analytical strong 
coupling behavior. The arrows indicate the expected asymptotic behavior for weak coupling which 
is proportional to vA, i.e. a constant independent of A in the plot. The solid lines show the actual 
analytic parameterizations in the weak coupling regime (cf. Eq. (|86|) ). 
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FIG. 6: Optimal wave functional parameter 5o(\,fj) as a function of A obtained from the simu- 
lation on a 16 3 lattice for three different values of fj 2 . The red shaded area corresponds to the 
phase transition region for all values of fj 2 . The dotted lines show the predicted analytical strong 
coupling behavior. The arrows indicate the expected asymptotic behavior for weak coupling which 
is proportional to \/A, i.e. a constant independent of A in the plot. The solid lines show the actual 
analytic parameterizations in the weak coupling regime (cf. Eq. (|86|i ). 
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0,0 0,2 0,4 0,6 0,8 1,0 

fj 

FIG. 7: Optimal wave functional parameter po(A, fj) as a function of fj obtained from the simulation 
on a 16 3 lattice for four different values of A. The expected A 1 / 2 71 (0) behavior for the equal time 
Hamiltonian with fj = 1 is scaled out. The solid lines show the analytical parameterizations. The 
dotted line corresponds to the "naive" analytical weak coupling prediction. 
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FIG. 8: Optimal wave functional parameter <5o(A, fj) as a function of fj obtained from the simulation 
on a 16 3 lattice for four different values of A. The expected A 1 / 2 71 (0) behavior for the equal time 
Hamiltonian with fj = 1 is scaled out. The solid lines show the analytical parameterizations. The 
dotted line corresponds to the "naive" analytical weak coupling prediction. 
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VII. CONCLUSIONS AND OUTLOOK 



Light cone coordinates are especially suited to parameterize high energy reactions for 
which perturbative QCD calculations have reached an unprecedented accuracy. In this paper 
we have addressed the question how to include non-perturbative features of QCD on the 
light cone. We propose to use lattice gauge theory formulated in exactly these coordinates. 

We start from the standard lattice action written in terms of near light cone coordinates 
such that the continuum action is recovered in leading order in the lattice spacing. The 
distance to the light cone is tuned by the adjustment of an external parameter r\. A transition 
to Euclidean time in this framework turns out to be problematic from a numerical point 
of view due to the fact that the Euclidean action remains complex which means that the 
integrand of the path integral cannot be interpreted as a probability measure anymore. 
Similar problems for QCD at finite baryonic density are generally referred to as the sign 
problem for which, up to the moment, no solution is known. In our case, this problem can 
be circumvented by applying the following strategy. We stay in Minkowski time, but switch 
to a Hamiltonian formulation of lattice gauge theory. Then the time evolution operator can 
be analytically continued to imaginary times and acts as a projector onto the exact ground 
state when it is applied to a trial state with a non-vanishing overlap with the exact ground 
state. Hence, instead of sampling the Euclidean path integral one manipulates a probability 
distribution for the product of the exact ground state wave functional and a guidance wave 
functional in a Quantum Diffusion Monte Carlo algorithm. For an improvement of the 
convergence of the diffusion Monte Carlo, the guidance wave functional should be sufficiently 
close to the exact ground state. The main goal of the present paper has been to develop a 
convenient and numerically realizable ground state projection operator and to propose such 
a guidance wave functional. 

We first work out the more obvious continuum formulation. The continuum near light 
cone Hamiltonian has an asymmetry in the longitudinal and transversal fields. The transver- 
sal fields are "enhanced" in the Hamiltonian in comparison to the longitudinal ones by a 
factor of 1/rf which is due to the underlying Lorentz transformation of the chromo mag- 
netic and chromo electric fields. Furthermore, the obtained near light cone Hamiltonian 
is similar to the classical Hamiltonian of a charged particle moving in an electromagnetic 
background field, which contains terms linear in the particle momentum. Such terms yield 
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complex branching ratios in a Quantum Monte Carlo algorithm which cannot be interpreted 
as probabilities and make it fail. However, this problem can be avoided. Linear terms in 
the QCD near light cone Hamilton operator can be compensated by the generator of lon- 
gitudinal translations. The QCD ground state is translation invariant, i.e. an eigenstate of 
the longitudinal translation operator with zero eigenvalue. Since the Hamiltonian commutes 
with the longitudinal momentum operator, the longitudinal momentum is not affected by 
time evolution. Therefore, one is able to construct an effective Hamiltonian feasible for a 
Quantum Diffusion Monte Carlo having the same ground state as the exact Hamiltonian by 
adding the longitudinal momentum operator to the exact Hamiltonian. 

Having checked feasibility in the continuum we derive the lattice Hamiltonian from the 
action via the transfer matrix method. We allow in general different lattice spacings in 
longitudinal and transversal directions. It is remarkable that the parameter 77 controlling 
the distance to the light cone multiplies the lattice anisotropy parameter £ which represents 
the ratio of the longitudinal lattice spacing and the transversal lattice spacing. Since these 
two parameters always appear together, there is no difference between the light cone limit 
and the anisotropic lattice limit. We can construct an effective Hamiltonian similar to 
the continuum case by adding the effective longitudinal momentum operator to the lattice 
Hamiltonian. 

We analytically compute the lattice ground state wave functional of the effective Hamilto- 
nian in the strong and weak coupling limit. In the strong coupling limit we obtain a product 
of single plaquette wave functionals similar to the equal time scenario. In the weak coupling 
limit, the solution is equivalent to the solution of the near light cone Hamiltonian with 
abelian fields, i.e. it is a multivariate Gaussian wave functional with a covariance matrix 
weighting correlations of field strengths at different spatial separations. 

Motivated by the strong and weak coupling solutions, we have constructed an effective 
ground state wave functional which smoothly interpolates between these two extreme results 
and which can be used as a guidance wave functional for a Quantum Diffusion Monte 
Carlo algorithm. It is a variational ansatz for the whole coupling range which contains 
a product of single plaquette wave functionals with two variational parameters. We have 
variationally optimized the parameters by minimizing the energy expectation value of the 
effective Hamiltonian with respect to the ground state wave functional. The effective ground 
state wave functional serves as a starting point for further qualitative explorations. It can 
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also be extrapolated to fj = and may be used to simulate correlation functions which 
appear in hadronic cross sections. 

The effective ground state wave functional can be improved by allowing also long range 
correlations in the wave functional. This is motivated by the observation that in the weak 
coupling limit, the covariance matrix elements of the analytical ground state wave functional 
which connect longitudinally separated spatial points become more and more important. 
An exponential ansatz then may contain plaquettes which are connected back and forth via 
long strings of gauge links. In the light cone limit the energetically most favorable string 
configurations are elongated along the minus direction. Such an ansatz may interpolate 
in the whole coupling range by allowing a covariance matrix with adjustable parameters. 



Numerical techniques [38| exist for a guided random walk in parameter space. So it may 
be possible to construct on the basis of an improved weak coupling solution a reasonable 
numerical procedure to obtain a good ground state for the effective lattice Hamiltonian. 

There have been strong advances in light cone physics recently in string and supersymmet- 
ric theory 39|, |40j. A careful study of near light cone theory in lattice QCD may supplement 
this successful work. 
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APPENDIX A: SOME USEFUL COMMUTATOR RELATIONS 

In this section we collect some useful formulae for the computation of matrix elements. 
First, we want to apply an arbitrary kinetic energy operator to the unperturbed strong 
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coupling ground state 



^0°^ ) (cf. Sec. IV Aj) multiplied by an arbitrary function of the links 



(0) 



(0) 



X)c,-[n?(j/),[n|(iO,/(OT)]] 



(0) 



(Al) 



The following double commutators are of special interest 



!W),Tr 



Re UkAx 



3,a,V 



[c k + ci) Tr 



Re( Uki(x 



(A2) 



and 



Tr 



Re C/ W (l) 



4 (c fc + q) 



Tr 



Ref £/ fc j(£ 



(A3) 



For the elementary plaquette, we have the following commutation relation 



II?($,Tr 



Ref t/fcz(x 



i Tr 
+i Tr 
— i Tr 
— i Tr 



y lit i { ?. r /,/(.r) 



2 

~2 



\m(ul{x)U kl {x)U k {x) 
Imf[//(f)C/ w (x)^(x) 



Imf C/ fc j(f 



(A4) 



In the following we assume an exponential ground state wave functional with exponent 
F({U}) where F({U}) is some arbitrary real valued functional of the links and 
unperturbed strong coupling ground state (cf. Sec. IV Al) 



^ 0) > is the 



v& ) = exp[F({[/})] 



(o) 



(0) 



exp[F({[/})] . 



(A5) 
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Then, the expectation value of the color trace of the momentum operator II" (y) squared 
with respect to the ground state Eq. (1A5I) is given by 



•y 



J,a,y 



(A6) 



The expectation value of the momentum operator II" (y) times an arbitrary functional 
G({U}) of the links is given by 



my) G({u}) 



*0 



-vl/ 



(0) 



(0) 



(0) 



(o) 



exp 



F({U}) m(y)G({U}) 



■ exp 



F({U}) 



(0) 



nj^.exp F({[/}) G({U}) 



exp 



(0) 



exp 



F{{U})\ G({U}) 
W($),e W \F({U}) 



(0) 



exp 



F({U}) G({U})Tl«i 



■ exp 



F({U}) 



(0) 



G({U}) W(y) 



(A7) 



APPENDIX B: DERIVATION OF THE TRANSFER MATRIX OPERATOR T 

In this appendix we construct the transfer matrix operator T propagating a spatial lattice 
configuration from one time slice to the next and the Hilbert space on which it acts. Op- 
erators are written explicitly in boldface. The Hilbert space on which T operates contains 
general states \^) which can be expanded in link states: 

|tf ) = J VU V(U) \U) . (Bl) 

The measure is VU in Eq. (1B1I) refers to the correspondent product of SU (2) Haar measures 

VU = Y[dUj(x) . (B2) 
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The inner product in this Hilbert space is given by 

= J t>u m\uy v{u) . (B3) 

We define the operator T such that its matrix elements in the link basis are given by the 
transfer matrix Eq. fl36|) 

(U(x'+) \TU(x + )^= T(x' + , x+) . (B4) 

The path integral for finite lattice of N T time slices with periodic boundary conditions can 
be written as the trace of the iV T -fold product of transfer matrices 

/il II dU 3 {x)e lS ^ = Ti[T N ^] . (B5) 

The transfer-matrix operator T is related to the Hamiltonian, the generator of time trans- 
lations 

T = e -ia+H ^ H = lim -— log(T) . (B6) 

a+^o ia + 

We define with the group elements gj(x) G SU(2) the following operators 

U 3 (x)\U) = U 3 (x)\U) Wj,x 
R(^(2)) \U) = \W) 

\W) = \..., 9j (x)U 3 (x),...) . (B7) 

Here all links in \W) coincide with the correspondent links in \U) except for the link Uj(x) 
which is left multiplied by gj(x). The operator H,(gj(x)) is similar to the translation operator 
in Quantum Mechanics. It is a unitary operator and satisfies the group representation 
property, i.e. 

R(^(£))R(^(f)) = R{gj(x) ■ g' 3 {x)) . (B8) 

The group elements gj(x) are parameterized by the exponential map which yields the Haar 
measure dgj(x) 

9j (x) = e^ )aa/2 , 7j a (x) G reals 
dg 3 {x) = j(7i(^))n^i( f ) ' ( B9 ) 
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The Jacobian J is equal to unity in a neighborhood of jj(x) = 0. Introducing the momentum 
operators II" (x) canonically conjugate to Uj(x) we have 



n«(x),uj,(f) 



e -i 7 ?0r)n?(a?) 



e i 7 ?(x)cr a /2 



(7, 



y Uj(x) S jtj > 8 3>3 > , 
-U](f) y 6 jd t 8 S)S > . 



(BIO) 



In contrast to the continuum commutation relations Eq. (1221) . the lattice momentum oper- 
ators canonically conjugate to Uj(x) do not commute. 



[[u^x),n b jt (x')],v k 



a o 
T' ~2~ 



1 e 



obc 



u fc (y) 



"j',fc ®x\y 

8x,x' 



= i e afec [lT](f), u 
Since this relation is true for arbitrary Ufc(y), we get 

[n^n^fO] = i s^ c u%x) 8 jif s. 



(Bll) 

(B12) 
(B13) 



Note that we have defined the translation operator on the group manifold R with an opposite 
sign inside of the exponential in comparison with 2]. Our definition yields the same com- 



mutation relations as 



34j which reproduce the continuum commutation relations Eq. (T22|) 



with the gauge field Aj(x) in the naive continuum limit. For simplicity we abandon to write 
quantum mechanical operators explicitly in boldface in the following. By using the group 
translation operators R we may write for the transfer matrix operator 



JJ J dg^(x) i2 (</_(£)} exp<j i— — ±- IV 
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1 [ / dg k (x) j 
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(B14) 



5( 2 a\ 



Tr 



1 - Re C/ 12 ff 



45 



It still has the right matrix elements Eq. (IB4I) . In order to arrive at Eq. (IB 141) one uses the 
fact that R(gj(x)) parameterizes the translation in group space from Uj(x, x + ) — > Uj(x, x' + ) 

g 3 (x) = Uj{x, x l+ )U}(x, x + ) . (B15) 

Now, one may perform the group integrations in Eq. (1B14I) explicitly. In the limit a + —>■ 0, 
the time evolution along one temporal step a + induces rotations 9j{x) which are of the order 
a + and are close to 1. This implies that the parameters jj(x) parameterizing these shifts are 
of the order a + as well. Therefore, it is convenient to make an expansion around 7" (it?) = 
up to order 0(a 2 + ). In this limit, the Jacobian is approximately equal to 1 and the integrals 
become Gaussian integrals which can be analytically computed. One obtains 

1 



lat 



lim 



9 



ia + 
2 1 



log(T) 



1 



k.a 1 



2 



cr, 



-j Im ( U -h(x) 



2 a 2 ^ " 



(x) 2 + A^Tr 
9 2 a± 



1 - Re U 12 (x 
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